Community detection in networks with positive and negative links 
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Detecting communities in complex networks accurately is a prime challenge, preceding further 
analyses of network characteristics and dynamics. Until now, community detection took into account 
only positively valued links, while many actual networks also feature negative links. We extend an 
existing Potts model to incorporate negative links as well, resulting in a method similar to the 
clustering of signed graphs, as dealt with in social balance theory, but more general. To illustrate 
our method, we applied it to a network of international alliances and disputes. Using data from 
1993-2001, it turns out that the world can be divided into six power blocs similar to Huntington's 
civilizations, with some notable exceptions. 

PACS numbers: 89.75.Hc, 89.65.-s 

Keywords: complex networks; community detection; modularity; negative links; social balance theory; inter- 
national relations 



I. INTRODUCTION 

Many complex phenomena can be represented as net- 
works, and subsequently be analyzed fruitfully [J, H, Q. 
One of the first targets of network analysis is the de- 
tection of communities on the basis of the links, i.e. the 
possibly valued, or weighted, arcs or edges that connect 
the nodes. Intuitively, an assignment of nodes to com- 
munities should be such that links within communities 
are relatively dense and between communities relatively 
sparse. This means we should compare actual densities 
to expected densities of links within and between com- 
munities. Furthermore, since nodes, for example humans 
or proteins, can be members of different communities at 
the same time, e.g. organizations or protein complexes, 
respectively, the assignment should allow for the possi- 
bility that communities overlap. 

In approaches to find appropriate community assign- 
ments, much progress has been made in recentyears 0, 
S @| by using a concept known as modularity fr\ . While 
current modularity approaches take for granted that links 
are positively valued, representing bonds or attraction, 
scientists in numerous fields grapple with networks that 
also have negative links that represent repel, conflict, 
or opposition, for example in neural networks, semantic 
webs, genetic regulatory networks, and last but certainly 
not least, in social networks. 

In this paper, we generalize an existing Potts model Q 
for positive links to incorporate negative links as well. 
We will follow the intuition that the assignment of nodes 
related by negative links should be done the opposite 
way of positive links, with negative links sparse within 
and more dense between communities, generalizing an old 
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idea from social balance theory [9| . Finally, we apply our 
approach to a network of conflicts and alliances between 
countries. 

Recently, it was shown that modularity mig ht miss 
small communities embedded in larger ones LlO], and IS 
less accurate if the actual communities are highly dif- 
ferent in size pT| . Our method has two balancing pa- 
rameters that address this problem to some extent [,12|. 
Yet community detection through modularity remains a 
global rather than a local approach. 



II. PROBLEM STATEMENT 

We consider a directed graph G with n nodes and m 
links, which can be easily generalized to weighted graphs. 
We denote the total number of positive links in G as to"'" 
and the number of negative links as to~, hence to = 
TO"*" + m^. We define the entries of the adjacency matrix 
of G as follows: if a positive link is present from node i to 
node j, Aij — 1, if a negative link is present, Aij = — 1, 
and Aij = otherwise. For a weighted graph the link 
values, or weights, are denoted by Wij. We separate the 
negative and positive links by setting A^- ~ Aij if Aij > 
and zero otherwise, and A'^- — —Aij if Aij < and zero 
otherwise, so A = — A^ . The positive and negative 
in- and outdegrees of i are defined as 

±fc°«* = ±fc™ = j:^a% (1) 

Our challenge is to assign each node i to one of c com- 
munities cFi € {l,...,c}. A complete configuration of 
community assignments is denoted by {cr}, which assigns 
each node i = 1, . . . , n to a community ai, . . . ,an- 
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III. SOCIAL BALANCE 

The challenge of community detection in networks with 
positive and negative links was first addressed by social 
balance theory, which has its origins in cognitive disso- 
nance theory from the 1950s. This theory is based on 
the notion that if two people are positively related, their 
attitudes toward a third person should match. For exam- 
ple, if Harry and Mary are positively related as friends, 
and both of them are related to John, they should both 
be related to him either positively or negatively. In either 
case their triad is said to be socially balanced. If Harry 
has a positive relationship with John while Mary is nega- 
tively related to John or vice versa, their triad is socially 
unbalanced. If all triads in a network are balanced, the 
network is said to be balanced. This definition was later 
generalized to cycles, a triad (a cycle of length 3) being 
a special case. 

The question whether a balanced network can be di- 
vided into separate parts arises naturally. The challenge 
is to define clusters of nodes such that there are only pos- 
itive link within clusters and negative links are between 
clusters. It was proven ^ that if a connected network is 
balanced, it can be split into two opposing clusters (and 
vice versa). 

However, there is an ambiguous case. If a triad has 
only negative relationships, it is neither balanced nor can 
it be split into two clusters. But it can obviously be split 
into three clusters. In order to accommodate for this 
possibility, the definition of balance was generalized to 
/c-balance [3, [3]. A network is fc-balanced if it can 
be divided into k different clusters, each cluster having 
only positive links within itself, and negative links with 
other clusters. It can be proven [l3| that a network is 
balanced if and only if it contains no cycles with exactly 
one negative link. The intuition is simple. Suppose there 
is a cycle viV2 ■ ■ ■ VkVi with one negative link, say between 
vi and Vk , and only positive links between the remainder 
nodes, then 1 and k are both positively and negatively 
connected, and the cycle is unbalanced. But if in this 
cycle there is also a negative link between i and j, and 
1 < i < j < k, we can split the cycle in two parts, one 
cluster from 1 to i and one from j to k. If there are 
more than two negative links, we can split up the cycle 
analogously into more clusters. 



IV. FRUSTRATION 

In reality, however, social networks are rarely, if ever, 
fully /c-balanced. The question then becomes whether 
we can still cluster nodes. Obviously, there are some 
links that make a network unbalanced. The number of 
such links can be expressed as an amount of frustration. 
Links that contribute to frustration are negative links 
within clusters and positive links between clusters. So 



the following expression should be minimized 

5^ AA-.<5(a., a,) + (1 - X)A±il - <5(a., a,)), (2) 

where S{(7i,aj) = 1 if at = (Jj and zero otherwise, and 
A is a parameter through which the contribution of the 
two types of links can be balanced. This is the approach 
taken in 

The objective, then, is to find a clustering {a} such 
that the frustration is minimized. Since the term (1 — 
A)A^ does not depend on the specific configuration and 
is therefore irrelevant for finding the minimum, we can 
simplify the above expression to 

5^(AA^.-(1-A)4)^(a„a,), (3) 

ij 

We can now see that only for A = 1/2 we retrieve A = 
^4+ — A^, up to a multiplicative constant of 2. Using any 
other value for A would change the minimum found, but 
changing A is the same as altering the (weights of the) 
original network. Setting A = 1/2 accordingly, we can 
simplify further, and now define frustration as 

F{{a})^-J2A^A<^^,a,). (4) 

However, frustration does not generalize to a network 
with only positive links. In that case, frustration groups 
together all nodes into one cluster. Even if there are 
some negative links, frustration will cluster together very 
sparsely connected nodes. It's therefore clear that this 
approach does not match with current methods of com- 
munity detection. Preferably, there should not be a dis- 
tinction between methods for positive and others for neg- 
ative links, but there should be one method for both. 

V. MODULARITY 

In approaches to find appropriate community assign- 
ments in networks with only positive links, much progress 
has been made recently [J S Q . The principal method 
for detecting communities is through modularity opti- 
mization, which boils down to clustering nodes based on 
the link densities within and between communities. The 
link densities should be high within communities and low 
between communities. 

The ordinary^ definition of modularity for directed 
graphs [l^ is 

QiW}) ^ —^niss - [rriss], (5) 



Keep in mind that we consider directed graphs. Therefore the 
sum of all degrees is m, not 2m, which is reflected in some minor 
changes to the original definitions [Tsll : see also |19|| . 
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a given configuration {cr}. We reward internal positive 
links by and penalize absent internal positive links by 



bij, which leads to 



c )k = -l 



FIG. 1: (Color online) Network illustrating the shortcoming of 
modularity when negative links are present. The dashed links 
have a weight of —1 and the others have a weight of 1. Using 



expected link values pij — 
equal the actual values, pij 
Q = for all configurations. 



k° k]"/m, the expected values 
— Aij. Hence the modularity 



where 



is the actual number of links within a com- 



munity s, J2 ^ijH'^ii 5)15(5, CTj), and [toss] is the expected 
number of such links, '^Pij5{ai, s)d{s,(jj), where pij is 
some expected value in a random null model. The ex- 
pected values are constrained by ^Pij = m, because the 
random null model should have the same number of links 
as the actual network. Taking degrees into account, a 
sensible expectation is pij = k°^'^kj"/m, which was used 
in the original definition [T9|."^ 

A straightforward generalization to weighted graphs 
is to set Aij — Wij and to take the degree measure as 
the sum of the link weights. For graphs with negative 
weights, however, a problem arises, illustrated in Fig. [1] 
The weighted degree of the three nodes a, b and c is 
ka = 1, kt = 1 and kc = —1. The total weighted degree 
is m = ^ Wij = 1 . The expected values pij 



equal the edge weights Wij. Hence Ai 



Pij 



k^^'^kf/m 
= for all 



links, and each possible community configuration results 
in a modularity Q = 0, while the appropriate configura- 
tion is clear from the figure: a and b belong to the same 
community, and c to another community [20| . Some 
adaptations of modularity are therefore required to de- 
tect communities in networks with (also) negative links. 



VI. POTTS MODEL AND EXTENSION 

To get rid of the shortcomings of frustration and mod- 
ularity, we will attempt to modify the latter. After all, 
for many networks with only positive links, modularity 
has shown to uncover valid community structures in mul- 
tiple areas . We will extend the original definition of 
modularity to allow for negative links, by phrasing our 
challenge as a Potts model [2'^. This model will also 
show up a close connection between social balance and 
modularity. 

We start out by treating the positive and the negative 
links separately. Mimicking the approach taken by Re- 
ichardt and Bornholdt [8| we first define a Hamiltonian 
for the positive part, which represents the "energy" of 



Setting ai 



1 — bii and 6, 



where pf^ repre- 



sents the expected (positive) link between i and j, allows 
us to simplify the above equation to^ 



n+i{a}) = -J2{A±-^+p+)6ia,,a,). 



(7) 



which is the Potts model analyzed by Reichardt and 
Bornholdt Q if only positive links are present. We de- 
fine the negative part analogously, but now we penalize 
internal negative links and reward absent negative inter- 
nal links. 



(8) 



The effect of the negative links on the energy of the entire 
configuration is opposite to the effect of the positive links. 
Combining the two Hamiltonians into one yields 

n{{a}) = (1 - x)n+i{a}) + xn-i{a}), (9) 

where A plays a similar role as in frustration, of balancing 
the effects of positive and negative links. As explained 
earlier, it makes sense to weigh the contributions of each 
part equally, thus A = 1/2. To illustrate, let us define a 
new matrix Bij = (1 — X)Afj — XA~j, and construct the 
Hamiltonian for this altered network by setting A — 1/2. 
Since the expected values for B are p'-'!^ = (1 — ^)ptj 
and p'^J = Xp~j, the Hamiltonian for B is equivalent to 
the one for A, up to a multiplicative constant of 2. So 
we may indeed set A = 1/2 and then simplify the above 
Hamiltonian (up to the multiplicative constant of 2) to 



H(W) 



E 



[A^J-h+pt,'J'p7J)]S{<J^,<T,), (10) 



which is the measure that we optimize to detect a com- 
munity structure in networks with both positive and neg- 
ative links. It can be easily seen that when the network 
is positive (and 7^ = 1) we obtain 



Q(W) = --H(W). 
m 



(11) 



So minimizing the Hamiltonian (|10p is the same as max- 
imizing modularity. In fact we just compare the origi- 
nal network to the appropriate negative link null model. 



Notice that if wc have a weighted network, aij = Wi. 
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which wasn't the case in the original modularity [18[ and 
in the Potts model [1]. 

The simplest version of the expected values, pfj, is ob- 
tained by just regarding the proportion of positive or of 
negative links in the network, p^- = /n(n — 1). If we 
want to take the degree distribution into account, then 
pfj = ^kf^* ^k'p/m^. The modularity given in [20] also 
defines this negative link null model appropriately, and 
is a special case of ours. 

When — 'y^ = 0, the Hamiltonian (fTU]) equals the 
frustration ^ of the network, and if the network is also 
balanced and complete (no missing edges), minimizing 
the Hamiltonian (|10[) yields the same result as minimizing 
the frustration (jl]). This can be pointed out by defining 
the probabilities by p^- = rri^ /n[n — 1), and by allowing 
the complete and balanced network to consist only of link 
values Aij G { — 1, !}• Then, as long as 7+™^ — j^m^ < 
n{n — l), the coupling Aij — {j'^p'lj ~7 Pij) between each 
positively associated pair of nodes is positive. Hence, the 
configuration produced by minimizing the Hamiltonian is 
the same as when minimizing the frustration. 

VII. THE GROUND STATE 

Finding the actual minimum of the Hamiltonian — the 
so called ground state — is NP hard [23], and therefore 
only heuristic methods can be applied. Our modular- 
ity (llip can be easily integrated with existing algorithms 
for modularity optimization, such as eigenvector [6j], ex- 
tremal optimization fast unfolding [23], or simulated 
annealing 0, [2^ . We opted for simulated annealing [2^ 
to minimize the Hamiltonian (fTU]) because it performs 
well in standard performance tests M. [27l. [28l [2^, al- 
though it's not the fastest algorithm (24l. l30jT^ Here we 
will give a short overview of how to adapt the simulated 
annealing approach to incorporate negative links. 

First, it's convenient to define the adhesion between 
community r and s, similar to Q, 

Urs ^ ("1+ - m,7j - (Jm+ ] - [m-^]), 

where m^^ — J2 ^fjH'^iT''")H'^jT ^) the actual number 
of arcs from r to s and [mfg] = '^J^pfjS{ai,r)d{aj,s) 
is the expected number of arcs from r to s. Hamilto- 
nian (jlOp can be rewritten accordingly, 

_ _ ^ (12) 

where r and s are communities in {a}. 

In order to minimize the Hamiltonian, we consider the 
effect of moving a single node from one community to 
another, expressed in terms of adhesion. More specifi- 
cally, moving node v from r to s results in the following 
change, 

/S.TL{(Jv'-r — > s) = (a„r + a™) - {a-vs + a™), (13) 



where the adhesion between node v and its comple- 

ment in community r. Let us write the mutual adhesion 
of a node v and a community r as Q!„(r) — a^r + oirv 
If, for the sake of argument, the mutual adhesion of v 
and s is larger than the mutual adhesion of v and r, 
A7i(tT„:r ^ s) = a„(r) — a„(s) decreases the Hamilto- 
nian. In other words, v has more positive, or less nega- 
tive, links than expected to s than to r, and moving v to 
s would improve the configuration. 

To each move we can assign a probability [26| 

exp(/3AH(CT„:r -> s)) 
exp(/3A7^((Ti,: r -> i)) 

where T is the temperature and (3 = 1/T. By slowly 
decreasing the temperature, the probability of moving to 
another state approaches the maximum possible, thereby 
forcing the system into its minimum energy, i.e. the 
ground state. Notice that in principle the probabilities 
are dependent on the total energy 7i((T„ : r s), but 
since 1-L{(Jv : r s) — T-L{{a}) + A'H{(Ty:r s) we can 
simplify to the equation stated above. 

The algorithm iterates randomly over the nodes a num- 
ber of times, after which the temperature is decreased to 
a lower temperature T' stepwise, and usually (although 
not necessarily) T' — 0.99T. The iterations and the low- 
ering of the temperature are continued until there are no 
further (significant) improvements. Any further changes 
would result in a higher energy, which we do not want, 
hence the resulting configuration of minimum energy is 
our solution {ct}. Herein, for any set of nodes u, its mu- 
tual adhesion to its own community s is stronger than 
to any other community r, a„(s) > au{r), which is clear 
when one looks at Eq. p^ . 

Furthermore, the cohesion, or self-adhesion, Uss of a 
community is always positive, ass > 0, and the mutual 
adhesion between two communities is always negative, 
ars+Qsr < 0. If the cohcsion wcrc negative, wc could then 
move a set of nodes to another community and thereby 
decrease the energy, which would contradict the fact that 
the system is in the ground state. 

In fact these last two inequalities can be rephrased, 
which yields some insight into the effect of the param- 
eters 7^. If we assume, for analytic purposes, that 
p,^ — p^ ~ m^/n{n — 1), the expected values become 
[m^s] = p^nrUs for r ^ s and [mfg] = p^ns{ns — 1), 
where Us is the number of nodes in community s. Writ- 
ing this out we arrive at 

iT^ts - '"^ss > 7+m+ - 7~m~ 
ns{ns - I) ~ n(n - 1) 

> (m+ - TO' ) -f (m+ - m" ) ^^^^ 

wherein the middle term is a sort of global density. Hence 
by changing 7^ we change the threshold for clustering 
nodes together versus keeping them apart. Either way, 
the density within a community is always higher than 
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the global density of the network, while the density be- 
tween communities is always lower than the global den- 
sity. Increasing 7+ raises the threshold for nodes to be 
clustered, and will (generally) result in smaller commu- 
nities detected, possibly embedded in larger and sparser 
communities. Increasing 7" has the opposite effect and 
lowers the threshold, and will (generally) result in a con- 
figuration of larger communities. 



VIII. APPLICATION 

To show how our method can be applied to an em- 
pirical network, we analyze international relations taken 
from the Correlates of War [3l|, [13] data set over the 
period 1993-2001, where military alliances can be repre- 
sented by positive links and disputes by negative links. 
The data set contains a wide variety of disputes, for ex- 
ample border tensions between Colombia and Venezuela, 
the deployment of Chinese submarines to Japanese is- 
lands, and Turkish groups entering Iraqi territory. Dis- 
putes were assigned hostility levels, from "no militarized 
action" to "interstate war," and we chose the mean level 
of hostility between two countries over the given time in- 
terval as the weight of their negative link. The alliances 
we coded one of three values, for (1) entente, (2) non- 
aggression pact, or (3) defense pact. The disputes w^j 
and alliances wf^ are both normalized to values in the 
interval wf^ € [0, 1] which then bear equal weight in the 
overall link value wa — wf, — w~-. The largest connected 
component consists of 161 nodes (countries) and 2517 
links (conflicts and alliances). 

The result of the analysis [Q = 0.561) is shown in 
Fig. [21 Countries of the same color (or pattern) belong 
to the same community, which in this context is more ap- 
propriately labeled a power bloc. How strongly a country 
belongs to its power bloc can be determined by the ad- 
hesion ay{s) a node has to its community. The power 
blocs can be identified as follows: (1) the West; (2) Latin 
America; (3) Muslim World; (4) Asia; (5) West Africa; 
and, (6) Central Africa. If we detect communities by us- 
ing only positive links, there is an agreement of about 
64% with the configuration in Fig. [21 while if using only 
negative links, there is an agreement of about 30%. 

Our result resembles the configuration depicted in 
Huntington's renowned book The Clash of Civiliza- 
tions 33], with a few notable exceptions. The West 
African power bloc is an additional insight that is ab- 
sent in Huntington's configuration. A major difference 
with Huntington is that China itself does not constitute 
a separate bloc, nor does Japan or India. Some other 
noteworthy differences are Pakistan and Iran which are 
grouped with the West, while South Korea and South 
Africa are grouped with the Muslim World. 

If we run the algorithm with 7+ = 0.1 and 7^ = 1, 
North America merges with Latin America, while Europe 
becomes an independent community, and North Africa 



and the Middle East align with Russia and China. When 
setting 7+ = 1 and 7" = 2, in contrast, former Soviet 
countries separate from Russia and form an independent 
community. Using a range of values for 7^, one can de- 
tect various layers in the community structure. 

Our configuration docs not imply that conflicts take 
place between power blocs only, as 24% of all conflicts 
actually take place within blocs. For example, Georgia 
and Russia had serious conflicts, and DR of Congo and 
Rwanda had theirs, but each of these pairs is grouped 
together nevertheless. In these cases, the alliances over- 
came the conflicts in the grouping, confirming that a con- 
figuration of international relations is more than the sum 
of bilateral links. 

Our political analysis here is limited, since we wish to 
demonstrate the method rather than present a complete 
coverage of international alliances and disputes. Other 
approaches that could be brought into play are the demo- 
cratic peace theory [sj, HH] , which predicts few conflicts 
between democratic countries but fails to predict that in 
actuality, most conflicts occur between democratic and 
non-democratic countries; the realist school [3^, which 
emphasizes geopolitical concerns; and, the trade-conflict 
theory [37| , which argues that (strong) trade relations di- 
minish the probability of a dispute, or lower its intensity. 
In sum, although Huntington's configuration of civiliza- 
tions was questioned ^38l . [3^ , it seems to be fairly robust 
and with some marked exceptions is confirmed by our 
analysis. 



IX. CONCLUSION 

We have extended the existing Potts model by adapt- 
ing the concept of modularity to detect communities in 
complex networks where both positive and negative links 
are present. This approach solves a long-standing prob- 
lem in the theory of social balance, namely the clustering 
of signed graphs. 

As a case in point, we have analyzed a social network 
of international disputes and alliances. Other applica- 
tions could be networks of references on the Web [i^l 
or in blogs SlJ. If in these data positive and negative 
references are distinguished, our method makes possible 
to detect not only thematic clusters, but also positional 
clusters with internal agreement and external disagree- 
ment. 

For network data, the model's parameters (7*) can be 
used to find smaller (sub) communities, although there 
is currently no theoretical guidance to choose parameter 
values 12]. Even if there were such guidance, the mod- 
ularity approach intrinsically aims at global rather than 
local optimization. Our implementation is based on sim- 
ulated annealing [2^ , which performs quite well with 
standard tests, although for very large networks, faster 
algorithms will be necessary [23 |. 

Whatever algorithms future researchers will use, or im- 
provements of the concept of modularity they will de- 




□ West H Asia 

Latin America s West Africa 

B IVIuslim World a Central Africa 

FIG. 2: (Color online) Map of the communities in the conflict and alliance network found using the algorithm described in the 
text (0 = 0.56f, 7=^ = 1). 



velop, being able to detect communities in networks with 
both positive and negative links is important in numer- 
ous fields of science, and a stepping stone toward further 
analyses of complex networks. 

Acknowledgments 

The authors like to thank anonymous referees and 
Jean-Charles Delvenne for their valuable comments and 



advice. VT acknowledges support from a grant "Actions 
de recherche concertes Large Graphs and Networks" of 
the "Communaute Frangaise de Belgique" and from the 
Belgian Network DYSCO (Dynamical Systems, Control, 
and Optimization), funded by the Interuniversity Attrac- 
tion Poles Programme, initiated by the Belgian State, 
Science Policy Office. 



[1] S. Boccaletti, V. Latora, Y. Moreno, M. Chavez, and 
D. U. Hwang, Phys. Rep. 424, 175 (2006). 

[2] S. N. Dorogovtsev, A. V. Goltsev, and J. F. F. Mendes, 
Rev. Mod. Phys. 80, 1275 (2008). 

[3] J. Bruggeman, Social networks: An introduction (Rout- 
ledge, London, 2008). 

[4] J. Duch and A. Arenas, Phys. Rev. E 72, 027104 (2005). 

[5] G. Palla, I. Derenyi, L Parkas, and T. Vicsek, Nature 
435, 814 (2005). 

[6] M. E. J. Newman, Phys. Rev. E 74, 036104 (2006). 

[7] M. Girvan and M. E. J. Newman, Proc. Natl. Acad. Sci. 
U.S.A. 99, 7821 (2002). 

[8] J. Reichardt and S. Bornholdt, Phys. Rev. E 74, 016110 
(2006). 



[9] F. Harary, Mich. Math J. 2, 143 (1953). 
[10] S. Fortunato and M. Barthelemy, Proc. Natl. Acad. Sci. 

U.S.A. 104, 36 (2007). 
[11] H. Du, D. R. White, Y. Ren, and S. Li, (to be published). 
[12] J. M. Kumpula, J. Saramaki, K. Kaski, and J. Kertesz, 

Eur. Phys. J. B 56, 41 (2007). 
[13] F. Heider, J. Psychol. 21, 107 (1946). 
[14] J. A. Davis, Hum. Relat. 20, 181 (1967). 
[15] D. Cartwright and F. Harary, Elemente der Mathematik 

23, 85 (1968). 

[16] P. Doreian and A. Mrvar, Soc. Networks 18, 149 (1996). 
[17] P. Jensen, Phys. Rev. E 74, 035101 (2006). 
[18] M. E. J. Newman and M. Girvan, Phys. Rev. E 69, 
026113 (2004). 



[19] E. A. Lcicht and M. E. J. Newman, Phys. Rev. Lett. 100, [31 

118703 (2008). 

[20] S. Gomez, P. Jensen, and A. Arenas, Phys. Rev. E 80, [32 
016114 (2009). 

[21] J. Reichardt, Structure in Complex Networks (Springer, [33 

New York, 2008). 
[22] F. Y. Wu, Rev. Mod. Phys. 54, 235 (1982). 
[23] U. Brandos, D. Dolhng, M. Gacrtlor, R. Gorkc, and [34 

M. Hocfcr, arXiv:physics/0608255v2 [physics. data-an] 

(2006). [35; 
[24] V. D. Blondel, J.-L. Guillaume, R. Lambiotte, and 

E. Lefebvre, J. Stat. Mech.: Theory Exp. 2008, P10008 [36 

(2008). 

[25] R. Guimera, S. Mossa, A. Turtschi, and L. A. Amaral, [37 
Proc. Natl. Acad. Sci. U.S.A 102, 7794 (2005). 

[26] S. Kirkpatrick, C. D. Gclatt, and M. P. Vecchi, Science [38 
220, 671 (1983). 

[27] L. Danon, D. A. Guilera, J. Duch, and A. Arenas, J. Stat. [39 

Mech.: Theory Exp. p. P09008 (2005). 
[28] R. Guimera and N. L. A. Amaral, Nature 433, 895 [40 

(2005). 

[29] A. Lancichinetti, S. Fortunato, and F. Radicchi, Phys. [41 

Rev. E 78, 046110 (2008). 
[30] M. E. J. Newman, Proc. Natl. Acad. Sci. U.S.A. 103, 

8577 (2006). 



F. Ghosn, G. Palmer, and S. Bremer, Conflict Manage- 
ment and Peace Science 21, 133 (2004). 

D. M. Gibler and M. Sarkees, J. Peace Res. 41, 211 
(2004). 

S. P. Huntington, The clash of civilizations and the re- 
making of world order (Simon & Schuster, New York, 
1996). 

P. R. Hcnsel, G. Goertz, and P. F. Diehl, J. Polit. 62, 

1173 (2000). 

D. Tocquevillc and A, Democracy in America (Perennial, 
New York, 2000(1835)). 

H. Kissinger, Diplomacy (Simon & Schuster, New York, 
1994). 

S. W. Polachek, J. Robst, and Y. C. Chang, J. Peace 
Res. 36, 405 (1999). 

E. A. Henderson and R. Tucker, Int. Stud. Q. 45, 317 
(2001). 

B. M. Russctt, J. R. Oneal, and M. Cox, J. Peace Res. 
37, 583 (2000). 

G. W. Flake, S. Lawrence, C. L. Giles, and F. M. Coetzee, 
Computer 35, 66 (2002). 

L. A. Adamic and N. Glance, in LinkKDD '05: Proceed- 
ings of the 3rd international workshop on Link discovery 
(ACM Press, New York, NY, USA, 2005), pp. 36-43. 



